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Abstract 

Two-temperature, two-mass quasi-equilibrium plasmas may occur in electron-ion plasmas, 
nuclear- matter, as well as in electron-hole condensed-matter systems. Dense two-temperature 
hydrogen plasmas straddle the difficult partially-degenerate regime of electron densities and tem- 
peratures which are important in astrophysics, in inertial-confinement fusion research, and other 
areas of warm dense matter physics. Results from Kohn-Sham calculations and QMC are used 
to benchmark the procedures used in classical molecular-dynamics simulations, HNC and CHNC 
methods to derive electron-electron and electron-proton pair-distribution functions. Then, nonequi- 
librium molecular dynamics for two-temperature, two-mass plasmas are used to obtain the pair 
distribution. Using these results, the correct HNC and CHNC procedures for the evaluation of 
pair-distribution functions in two-temperature two-mass two-component charged fluids are estab- 
lished. Results for a mass ratio of 1:5, typical of electron-hole fluids, as well as for compressed 
hydrogen are presented. 
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I. INTRODUCTION 



The study of hot strongly-coupled dense charged fluids is a difficult task, especially near 
the regime of or molecular and atomic species [l], or excitons in electron-hole plasmas. The 
system is better understood for fully ionized systems, such as hydrogen, in the form of 
free electrons and protons, and fully- ionized electron-hole condensates. In fact, considerable 
headway has been made using methods based on density-functional theory (DFT), even for 
plasmas with multiple states of ionization. DFT methods have been used with molecular- 
dynamics based approaches^ B, Q, as wel. as within mnhi-eomponent integral-eonation 
approaches [5|, [6[ . Equilibrium properties of plasmas, as well as their linear transport prop- 
erties, have been successfully studied in these papers, and excellent agreement between the 
molecular-dynamics based DFT and integral-equation based DFT has been found Q]. 

On the other hand, laser-produced plasmas are initially formed as two-temperature plas- 
mas, where the electrons have absorbed the laser energy and have self-equilibrated to some 
"electron temperature" T e , while the ions remain cool, at some temperature Tj, with Tj < T e . 
The opposite situation arises in shock-wave generated plasmas, where the ions absorb the 
shock energy and T; > T e . Such two-temperature plasmas also occur in astrophysical set- 
tings, affecting the time of termination of synthesis of light-nuclei to occur at different stages 
of cooling of the electrons^], and influencing the Coulomb nuclear-tunneling rates {{J. The 
possibility of such well-defined two-temperature plasmas is largely a result of the extreme 
mass ratio mj/m e > 1836 between ions and electrons. Similar, but less well defined situa- 
tions can arise in electron-hole plasmas, where the masses are of the same order of magnitude 
(e.g, the electron and hole masses in GaAs are 0.067m e and 0.34m e respectively, with an 
electron-hole mass ratio of ~ 5). GaAs is a direct bandgap material, and electron-hole 
plasmas are more easily studied in indirect-gap systems like Si where the density-of states 
mass ratio is ~ 3. The simulation of such systems at two temperatures, using quantum 
Monte-Carlo methods is at present unavailable, even in regimes of densities and tempera- 
tures where bound states (or exciton formation in electron-hole systems) do not exist. Thus 
it is natural to look for analytical methods based on integral-equation techniques which 
are computationally simple and physically insightful. However, although T e , T, define the 
temperatures of each subsystem and the pair-distribution functions (PDFs) g ee and ga, the 
'temperature' T e j entering into the cross-correlations g e i, as well as the the effects of electron 
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spin, exchange etc., relevant to two-temperature systems need to be clarified. In this context 



we use T e . 



T T- 

-*- Pi J- 1.0. 



Ti, and T e j to refer to the electron-, ion-, and electron-ion temper- 



atures as they enter independently into the Ornstein-Zernike (OZ) and hypernetted chain 
(HNC) equations. In fact, some authors 10J have proposed to modify the well-established 
OZ equations in dealing with two-temperature (2T) two-mass (2M) systems. 

The objective of this paper is to study such 2T-2M plasmas u sing results from molecular- 
dynamics (MD)|ll|, HNc(l2j, classical-map HNC (CHNC) [IJ, bij , quantum Monte Carlo 
(QMC) ll(] and Kohn-Sham (KS) 6] methods to establish the proper implementation of 
quantum effects and 2T, 2M situations in simulation studies. One of our main interests 
would be uniform hydrogenic plasmas free of bound states, in the regime of warm-dense 
matter. 



II. THEORETICAL METHODS 



A system of classical particles, e.g., hard-spheres or Lennard- Jones fluids, or classical 
ions in a uniform neutralizing background, can be studied completely using the method of 
molecular dynamics (MD) where the classical equations of motion are integrated numerically, 
for a sufficiently large number of particles contained in a simulation box. It has been 
found that the particle distribution functions, e.g., the pair-distribution function (PDF) 
9ij{ r ) (where i,j specify the species, spin etc.), obtained from MD simulations for charged 
classical ions can also be accurately reproduced via suitable integral-equations which are 
computationally very economical and efficient. The pair-potentials, quantum corrections 
etc., needed to simulate systems with ions and electrons, or purely electron systems (with 
a uniform neutralizing background) will be discussed in this section. First, we compare the 
usual HNC approximation with the CHNC method. Then, these results are compared with a 
full quantum (Kohn-Sham) calculation; these results are then used to determine the effective 
diffractive interaction used in the CHNC method. Finally, we discuss how the CHNC can 
be extended to "classical map molecular dynamics" (CMMD). 
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A. HNC and CHNC methods 



The HNC equation and its straight-forward generalizations, coupled with the OZ equation 
have lead to very accurate results for classical charged-particle interactions. The exact 
equations for the PDFs are of the form: 

g^(r) = e~^ cf ^ i3 ^ +hi3 ^~ Ci3 ^ +Bi3 ^ (1) 

Here <pij{r) is the pair potential between the species If the bridge function Bij(r) is 
set to zero we have the HNC approximation. Then, given the temperature T = 1/13, the 
particle densities n^, and the pair-potentials the pair-correlation function /iy(r) and 

the direct correlation function Cy (r) can be self-consistently obtained via the HNC and OZ 
equations, which have the form 

hij(r) = cy(r) + ^n s dr'h is (\r - r'\)c sj {r'). (2) 

s J 

This method already fails for strictly attractive potentials. Thus, in simulations of electron- 
proton systems, e p(?") has to be replaced by effective-potentials which attempt to incor- 



porate quantum diffraction effects 



151 ]. Even with purely repulsive potentials, classical 



simulations fail to incorporate Fermi or Bose statistics which begin to manifest as the den- 
sity is increased and the temperature is lowered. A well established means of incorporating 
quantum effects is to derive integral equations from correlated-determinantal wavef unctions, 
as done in the Feenberg approach {if]. u\ . The resulting integral equations are very daunt- 
ing, and in fact, not easy to use. Quantum Monte-Carlo (QMC) itself may be considered as 
an adaptation of the Feenberg-functional to generate a statistical measure for the stochastic 
algorithms used in MD. An alternative approach, using Feynman paths instead of classical 
trajectories, provides another class of simulation techniques. However, these quantum simu- 
lation methods become computationally extremely heavy. Such methods are best suited for 
the establishment of bench-mark results, and for the "calibration" of other methods which 
contain approximation schemes. In fact, the QMC techniques have been most useful in 
providing the "exchange-correlation" potentials V xc (r) needed in the Kohn-Sham density- 
functional theory (DFT) equations. While DFT is itself exact, one has to use results from 
QMC and such microscopic methods to model the unknown V xc (r). Given the V xc , the 
inhomogeneous density distribution around a given particle can be calculated, and the pair- 
distribution is deduced from it. 
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The method followed here is to exploit the well established Kohn-Sham equations as 
the reference calculation, and determine the effective potentials to be used in the classical 
simulations of 2T and 2M systems. To this end we present comparisons of Kohn-Sham 
calculations of g ep {r) for H-plasmas with available QMC results to mutually validate these 
methods. Another theoretical tool we use is the "classical-map HNC", i.e., CHNC equa- 
tions which incorporate quantum effects including Fermion statistics via effective potentials 
and effective temperatures. The CHNC has been extensively tested via comparisons with 
QMC results in 2-D and 3-D electron systems, and shown to provide excellent agreement, 



even at the extreme quantum limit of zero temperature 13|, ll8|] . CHNC uses a "quantum 



temperature" T q , which depends on the Fermion density. If the physical temperature of the 



quantum fluid is T, the distribution functions are obtained 
temperature T c f such that: 



13| from a classical fluid at the 



T cf = yjT 2 + T*. (3) 
The temperature T is defined to be such that the classical Coulomb fluid has the same 

n 

exchange-correlation energy as the quantum fluid [13J. DFT assures us that only the true 
density distribution possesses the true exchange-correlation energy. Thus, the charge distri- 
butions, i.e., the PDFs obtained from CHNC are found to be in excellent agreement with 
those from Monte-Carlo simulations of 2D and 3D systems. This agreement is obtained 
by including the exchange-hole of parallel-spin electrons as an effective potential, called in 
CHNC the Pauli exclusion potential Pij{r). Clearly, this is zero ifi^j. For i = j = || spins, 
the potential Pij{r) is such that the non- interacting PDFs, i.e., g^Ar) are correctly recovered 
from the integral equations [19]. Thus, using atomic units where |e| = K — m e = 1, the 
effective pair potentials <pij{r) are of the form: 

Mr) = Prjir)+V^{r) (4) 
V^r) = ZlZj (l- e - k ^)/r (5) 

Here z e = — 1, z p = 1 and kij is a cut-off 'momentum' defining a diffraction correction 
allowing for quantum effects. In the simplest formulation kij is the thermal de Broglie 
momentum given by: 

% = kff = (27rm ij T lj ) 1 / 2 (6) 
rriij = miiTij / (rrii + rrij) (7) 



The temperatures entering into the HNC equations for the gij(r) are given by: 

Tij/rriij = Ti/ m, + Tj/ rnj. (8) 

The large mass of the proton ensures that the diffraction correction, as well as the T q , is 
negligible for the proton-proton scattering process. A more complete approach to deter- 
mining kij is to solve the corresponding Kohn-Sham equation for the two particles in the 
Kohn-Sham potential of the medium, and matching the so that the quantum and classical 
values of the PDF at contact are in agreement. 

The Pauli exclusion potential is usually determined by inverting the HNC equations 
applied to the exactly known non-interacting quantum PDFs g%{r) of the uniform electron 
fluid. Here i runs through e |,e J., i.e., a spin-resolved, two-component electron system is 
used. If i 7^ j, Pij = 0. In the absence of strong magnetic fields, the spin- resolution is not 
needed in warm dense systems. Treating the electrons as a spin-averaged one-component 
subsystem simplifies the simulations. Due to the non-linearity of the inversion of the HNC 
relations given in Eqs. [9J the Pauli exclusion potential for paramagnetic electrons is not a 
simple average of the Pauli exclusion potentials of the spin-resolved cases. The corresponding 
Pauli potential P e (r) has to be extracted directly from the averaged g® e {r). Thus we have, 
for the spin-resolved and -unresolved cases: 

(3P u (r) = -H/S(r)]+JY2(r), (9) 
(3P e (r) = - H{gl(r) + l}/2] + iV°/2 (10) 
*2(r) = *&(r) - 4(r). (11) 

A bridge term (r) is used to correct the HNC for multi-particle effects poorly rendere d by 



HNC. Such bridge corrections are found to be very significant in the 2D electron fluid 18( , 
but not for 3D electrons at densities and temperatures considered in this study. 

The CHNC differs from HNC in the use of the Pauli potentials and the quantum tem- 
perature T q when treating quantum fluids. Also, the pair-potentials used in HNC have 
been constructed to agree with KS-charge profiles (see below). Hence any insights ob- 
tained for the two-temperature two-mass HNC can be easily transfered to the CHNC. The 
two-t emp erature electron-ion plasma was discussed in a formal analysis by Boercker and 
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201 ]. using a product form for the partition function. However, no comparisons of their 



results with actual simulations are available. The more general two-temperature two-mass 



HNC type equations have been discussed, most recently by Seuferling et al. n|. Using an 
analysis based on the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy as well 
as some factorization assumptions, the authors of Ref. [h]] have proposed modified OZ type 
equations for 2T-2M plasmas. While their formulae reduce to the usual OZ equations, viz., 
Eq. [21 for the m a » m b , T a = T b , the case m a = m b , T a = T b is not correctly recovered. 
The results presented in our study imply that the usual OZ equations hold in all cases, as 
long as the correct mass-dependent Tij, Eq. [8] is used in 2T-2M systems. 

B. The Kohn-Sham reference calculation 



Kohn-Sham theory at finite temperatures 6] states that the true one-particle density 
distribution of the system subject to an external potential is such that the free energy of the 
system is minimized. This theorem holds rigorously for a system in equilibrium and we use 
it to derive distribution functions by considering the inhomogeneous electron distribution 
around a proton in the plasma. Let n(r) and p(r) be the electron and proton charge densities 
around the proton at the origin. These tend to the average densities n = p far away(r — ► oo) 
from the proton at the center; then, 

g ep (r) =n{r)/n. (12) 

Instead of using a two-component DFT procedure, we make the further approximation, 
well established in practice, where the proton subsystem is replaced by a uniform positive 
background with a cavity, viz., a Wigner-Seitz sphere of radius r s = [3/ (47m] 1 / 3 . The 
positive charge scooped out to form the cavity is placed as a point charge at the origin, and 
forms the central proton. The finite-temperature Kohn-Sham equation is a consequence of 
the Euler equation for the stationary property of the free energy under functional derivation 
with respect to the electron-density distribution, viz., 

« = a (13) 

on(r) 

A standard Kohn-Sham type analysis now leads to the equation: 

[-V 2 /2 + Z/r - V ks (r)}Mr) = ( 14 ) 

where 

V ks (r) = V p (r,n(r)) +V xc (r,n(r),T e ) 



n(r) = Y,\Ur)\ 2 f{e v /T e ). 

V 

Here 

Vp(r,n(r)) = J dr'n[r')/\r — r'| 

is the Poisson potential of the electron distribution n(r). This distribution is evaluated 
self-consistently from the Kohn-Sham wavefunctions ip u (r), v — n,l, m, energy e u , with the 
occupation factor given by the Fermi function f(e/T e ). The potential due to the proton 
at the origin is Z/r, with Z = 1, and V xc (r,n(r),T e ) is the finite-temperature Kohn-Sham 
exchange- correlation potential U{ which depends self-consistently on the charge profile n(r). 
This is evaluated using the local-density approximation (LDA), unlike in CHNC where a 
fully non-local V xc (r) is evaluated via a coupling-constant integration over the g ee {r). The 
Kohn-Sham procedure uses only n{r) = ng ep {r), and does not provide a g ee (r). Since this 
problem contains only one proton, there is no proton temperature in the theory. However, 
due to the large mass of the proton, and due to the exclusion of other protons by the 
central proton (modeled by the Wigner-Seitz cavity), the value of g e p[f) at r — > given by 
the Kohn-Sham calculation is expected to be a valid estimate for the full electron-proton 
plasma. In fact, in the two-temperature electron-proton plasma, T ep of Eq. [8] reduces to 
T ee , as in the Kohn-Sham calculation, since m p » m e . That this one-proton Kohn-Sham 
calculation correctly reproduces the g ep {f) of the plasma, even at very low temperatures, is 
seen from the comparisons given in Fig. (TJ where the path-integral Monte Carlo (PIMC) 
PDFs for hydrogen from the work of Militzer and Ceperley 2l(] have been used. 



C. The effective electron-proton interaction for classical simulations. 

In a classical simulation of an electron-proton plasma, or in a CHNC calculation, the 
Coulomb interaction V ep (r) appears. This is the attractive classical potential associated 
with the quantum-mechanical operator —l/r ep . It is expected to have the form: 

/r (15) 
Kp = k^p f ep (16) 
k% = [2vrT ep m ep ] 1/2 (17) 

The thermal de Broglie momentum k dB provides a first approximation to k ep . But we 
choose k ep such that the g ep {r — ► 0) generated by the classical procedure (e.g., MD or 
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Vepir) 




FIG. 1: (Online color)Panels (a),(b) present a comparison of the DFT, QMC, CHNC, MD, and 
HNC calculations of the electron-proton PDF. The last three use the effective potentials of Eq. [T5l 
A quantum temperature T q , a Pauli potential and the effctive potentials with f ep different from 
unity are used in CHNC. The lower panels compare the spin-resolved electron-electron PDfs in the 
H-plasma, obtained from QMC and CHNC. P-Dw(DFT) and CHNC calculations use the formula- 
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tions of Dharma-wardana and Perrot 

|a 1131 • Th e QMC (PIMC) is from Militzer and Ceperley 
while the H-M(HNC) follow the MD calculations of Hanson and MacDonald[23] using the potentials 
of Eq. M 



CHNC) agrees with the g ep (r — » 0) obtained from the Kohn-Sham calculation at the given 
r s and T e . It turns out that the correction factor f ep is quite close to unity for sufficiently 
high temperatures. Even at T=10.77 eV, r s = 1, i.e., T/E F — 0.215, / e =0.922 and we see 
from Fig. [1] that the agreement between QMC, DFT, and CHNC is quite good. We have 
determined the value of f ep as a function of T e , r s by matching the CHNC calculation and the 
Kohn-Sham calculation. In effect, f ep is similar to a pseuodopotential or form factor for the 
electron-proton interaction. When bound states begin to be formed (r s > 1.8), the form of 
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f ep becomes more critical, but this problem does not arise within the densities studied here. 
However, classically, for attractive potentials, dynamical instabilities could occur at any r s , 
T and these have to be controlled using close-approach cutoffs on the potentials, as well 
as controls on the velocity distribution functions, to maintain the meaning of "subsystem 
temperatures" which are set to T e and Tj. Investigation of such instabilities where the 
velocity distributions do not conform to the two-temperature model is outside the scope of 
this study. 

The electron-electron interaction used in CHNC, and MD simulations, is also a diffraction- 
corrected Coulomb potential, V ee {r), with k ee being (27rm ee T c j) 1//2 , with m ee = m e /2 and 
requiring no additional correction factors. This diffraction correction can be derived from 



the Schrodinger equation describing electron-electron scattering 



2J. 



D. Classical-map Molecular dynamics 

The HNC method using all three items: (i) diffraction-corrected effective potentials, (ii) 
the Pauli exclusion potential and (iii) the quantum temperature T q , is the CHNC scheme. If 
the same three items were included in classical molecular dynamics simulations, we have a 
classical-map-MD scheme (CMMD). The CMMD is superior to CHNC since it automatically 
includes any bridge corrections, etc. that are not in the HNC scheme. In CMMD the electron 
temperature would be T c f, Eq.[3], as in CHNC. However, since bridge corrections are expected 
to be negligible in the regime of densities and temperatures considered here, we do not carry 
out CMMD simulations. 



III. RESULTS 



We first provide comparisons between simple classical MD simulations and HNC calcu- 
lations of PDFs of 2T-2M systems using the simplest diffraction corrected pair potentials, 
given by: 

<f>% = Zl z 3 {l-e- k ^)/r (18) 
kf = [2™^,] 1/2 (19) 

The MD simulations only need the individual subsystem temperatures Tu, Tjj, and no cross- 
species temperature T^, % ^ j is needed. This is achieved by employing two velocity-scaling 
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thermostats that adjust the electron and ion velocity distributions to have the desired mean 
values. In contrast, the HNC needs a specification for T^. Seuferling et al.flO] have suggested 
that the OZ equations also need to be modified. These issues can be tested by comparison 
with the MD results. 



A. Two-mass two-temperature systems. 

Systems where the two masses m a and mb of the two components a, b are equal cannot 
produce two-temperature quasi-equilibrium systems unless V a b is, for some reason, extremely 
different from V&& and V aa . Thus two-temperature plasmas may exist for significant times, 
even when the mass-ratio is of the order of 3-10, as in some solid state electron-hole plasmas 
where band-structure effects associated with the existence of indirect gaps introduce restric- 
tions on electron-hole recombination. Here we present HNC calculations of the PDFs of 
plasmas with mb/m a = 5, and compare them with MD simulations, to establish the correct 
implementation of HNC and CHNC procedures. 

In Fig. [2] we show the PDFs calculated for a two-component system with a mass ratio of 
5, using both the HNC with the standard OZ relations and MD. The MD simulations used 
300 particles, 40,000 equilibration steps, a time step of 0.02 of the inverse electron-plasma 
frequency, and data was then accumulated over 120,000 steps using the two thermostats 
described above. In the HNC calculation, the pair-potentials are given by Eq. [TBI an d the 
cross-species temperature is as in Eq. [Bl This simple HNC-OZ procedure is in very good 
agreement with MD, both for the equilibrium and non-equilibrium (two-temperature) cases, 
and we conclude that the additional procedures proposed by Seuferling et al.[lo| in their 
Eq. (37) are not needed. That is, our results show that a modified OZ equation is not 
necessary. The comparison between the HNC and the MD establishes the correctness of 
the basic HNC procedures even in the quasi-equilibrium case where the formal derivation of 
the HNC equations becomes an open question. However, once the correct HNC procedure 
is established, the calculations for the quantum two-temperature two-mass system can be 
carried out using the CHNC, with the same temperature assignments extended to include 
the quantum temperatures T q , and the Pauli potentials. 
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FIG. 2: (Online color)The upper panel shows the PDFs for a two-component (a, b) system at a 
temperature T a = Tf, = 100 eV, r s = 1, the masses m a , nib being 1 and 5. The component a 
electron like, while b is hole- like, i.e., positively charged. Lower panel shows the two-temperature 
two-mass case, with T& lowered to 30 eV. In both cases the PDFs g%j{r) calculated using the 
standard HNC and the OZ equation, Eq. [2] agree well with the MD results using the same input 
potentials as in HNC. 



B. Electron-proton systems in thermal equilibrium 

The simple diffraction corrected potentials, Eqs. [T8lwere used by Hanson and MacDonald 
in H-plasma simulations 231] . Their PDFs can also be generated using the simple HNC 
equations if the above 4>ij(r) are used. Hence, in Fig. [T]we have labeled the corresponding 
g ep (r) as H-M(HNC,MD). The PDFs obtained from the DFT calculation, using Eq. El as 
implemented in the codes by Perrot and Dharma-wardana 22], as well as the PIMC results 
of Militzer and Ceperley are also shown, to establish that these two first-principles methods 
are in excellent agreement. Here we note that the CHNC results for g ep and also the spin- 
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resolved g ee are in excellent agreement with the QMC PDFs. To obtain this agreement, the 
CHNC uses the slightly modified diffraction parameter k ep = f ep k^ p with / ep =0.922 and 
0.965 at r s =l and 1.5 respectively, as obtained by matching the contact value of the CHNC 
g ep to the value from the KS calculation. The MD-HNC using the Hanson-MacDonald 
approach leads to a large value of g ep at r — > 0, while the g ee (not shown in the figure) are 
in strong disagreement. The agreement between QMC and CHNC shown in Fig. [1] holds 
even better at higher temperatures, and this justifies our use of the CHNC and Kohn-Sham 
results as the reference calculations when QMC results are not available. 

C. Two-temperature electron-proton systems 

In this sub-section we compare classical two-temperature H-plasmas and show that the 
temperature T ep that appears in the cross-species HNC equation is indeed the electron 
temperature T ee , as in Eq. El for the limit m p » m e . Thus, we use the same T ep in the 
CHNC, to include the quantum corrections and compute gij{r). In Fig. [3]we show the cross- 
species (electron-proton) PDF for a hydrogen plasma with T e = 100 eV, r s = 1, for the four 
cases: Tj = 100, 60, 30 and 10. From panel (a) we see that the classical procedures (HNC 
and MD) using the simplest set of (f>ij{r), Eq. HH overestimate the g ep in comparison to the 
Kohn-Sham (DFT) estimate. In panels (b-d) we have two-temperature plasmas, and the 
MD calculation (which needs only T e , Tj) is closely reproduced by the HNC if T ep is set to 
T e . In panel (b) we show that the choice T ep = (T e + T p )/2 in the HNC is clearly inapplicable 
if the system is entirely specified by T e , T p , and r s . As seen from Fig. [3], quantum effects 
may significantly modify the PDFs even when the electrons are at 100 eV. 

Thus, in Fig. HI we present CHNC calculations for a two-temperature plasma with T e = 
100, Tj = 30 at the density r s = 1. The top panel shows that the proton-proton PDF 
calculated from the quantum procedure (using CHNC) is more strongly coupled than in the 
classical (using HNC) g pp . The stronger e-p interaction in the classical system, as shown 
in the enhanced g ep , leads to greater screening, weakening the ion-ion interaction. The 
lower panel shows the spin-resolved e-e PDFs, labeled g uu , g u d obtained from CHNC, and 
the classical g ee obtained from HNC. The CHNC correctly incorporates the exclusion effects 
via the Pauli potential, Eq. 
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FIG. 3: (Online color) Dense hydrogen: Panel (a) compares HNC and MD g e p(r) using the simplest 
set of classical potentials (Eq. ITS]) . The DFT PDF shows that the classical potentials are an 
overestimate. In panels (b-d) we use the same classical potentials to establish that the temperature 
T ep needed in the HNC is indeed T e if HNC and MD are to agree for two-temperature electron- 
proton systems. 

D. The electron-proton PDF and Pauli exclusion effects 



The electron-proton pair distribution function is mainly determined by the e-p interaction 
which is spin- independent. However, once an electron is correlated with a proton, the 
correlation of that electron with other electrons would be affected by Pauli exclusion effects 
associated with the electron spin. In the CHNC and CMMD schemes, the effect of the 
Pauli principle are incorporated as a potential, Eq. O between parallel-spin electrons. This 
potential is not used in MD and in the pure HNC scheme. Hence the g ee (r) obtained from 
HNC, shown in the lower panel of Fig. IH(b) is identical for parallel and anti-parallel PDFs. 
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FIG. 4: (Online color)The upper panel shows the e-p and p-p PDFs for an electron-proton plasma 
with T e = lOOeV, T p = 30eV, r s = 1, calculated using HNC and CHNC. The lower panel shows the 
e-e PDFs, where the HNC does not incorporate the effects of the exclusion principle. The CHNC 
g uu and g u d refer to spin parallel and antiparallel PDFs respectively. 

However, Fig. H^a) shows that the g ep obtained by the full CHNC, inclusive of the Pauli 
potential, f ep , and T q features is quite close to the pure-HNC result where f ep = 1 in the 
diffraction potentials. At r s = 1, T q /Ep = 0.768, and hence, when T e = 100 eV, i.e., 
T e /Ep = 1.9956, then T q itself is substantial. Thus the larger value of g ep (r) at r — > 
found in the HNC and MD is not due to the Pauli exclusion effects, but due mainly to two 
reasons: (i)the overestimate contained in the zeroth set of effective potentials where f ep = 1, 
and (ii)the use of the physical temperature T e as the effective temperature of the classical 
electron fluid, while T c f > T e is used in the CHNC. To check these, we have run CHNC 
calculations where (i) the Pauli potential was switched off while the T q , f ep were included; 
(ii) only the Pauli and f ep were included; (iii) only the Pauli and T q were included; and so 
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forth. Such "numerical experiments" enable us to conclude that the Pauli exclusion effect 
is of relatively low importance for the g e p{f) when T e is 100 eV and r s = 1. 



IV. CONCLUDING DISCUSSION 

The simplest classical rendering of quantum plasmas, based on the use of diffraction 
corrected potentials (Eq. [IS]) was used with HNC calculations and MD simulations to resolve 
the ambiguities and difficulties in handling the two-temperature, two-mass system. We 
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conclude that the modifications to the OZ equations proposed by Seuferling et al. 
not needed. The classical mapping of quantum systems to the HNC equations, as used 
in the CHNC was confirmed by comparisons with Kohn-Sham DFT calculations as well 
as with available PIMC results for compressed hydrogen plasmas at finite temperatures. 
We conclude that the HNC and CHNC, together with the standard OZ equations provide 
excellent, accurate and simple analytical tools for the investigation of many-particle quasi- 
equilibrium systems for which direct quantum simulations continue to remain too prohibitive 
or unfeasible. 
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